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Abstract. Photoproduction of 77 mesons off 12 C, 40 Ca, 93 Nb, and nat Pb nuclei has been measured with a 
tagged photon beam with energies between 0.6 and 2.2 GeV. The experiment was performed at the Bonn 
ELSA accelerator with the combined setup of the Crystal Barrel and TAPS calorimeters. It aimed at the 
in-medium properties of the Sn(1535) nucleon resonance and the study of the absorption properties of 
nuclear matter for r\ mesons. Careful consideration was given to contributions from r/n final states and 
secondary production mechanisms of 77-mesons e.g. from inelastic nN reactions of intermediate pions. The 
analysis of the mass number scaling shows that the nuclear absorption cross section <tjv^ for r\ mesons 
is constant over a wide range of the t) momentum. The comparison of the excitation functions to data 
off the deuteron and to calculations in the framework of a BUU-model show no unexplained in-medium 
modifications of the Sn(1535). 

PACS. 13.60.Le Meson production - 14.20.Gk Baryon resonances with S=0 - 25.20.Lj Photoproduction 
reactions 



1 Introduction 

The study of possible in- medium modifications of the prop- 
erties of hadrons is a challenge for both theory and exper- 
iment. In contrast to any other composite system, most of 
the mass of hadrons is generated by dynamical effects from 
the interaction of the quarks. An important role is played 
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by the spontaneous breaking of chiral symmetry, the fun- 
damental symmetry of QCD. The symmetry breaking is 
reflected in a non-zero expectation value of scalar qq pairs 
in the vacuum, the chiral condensate. However, model cal- 
culations (see e.g. Ref. pQ) indicate a temperature and 
density dependence of the condensate which is connected 
to a partial restoration of chiral symmetry. In this way, 
hadron in-medium properties are closely connected to the 
non-perturbative aspects of low-energy QCD. While a di- 
rect relation between the quark condensate and the in- 
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medium masses and widths of hadrons is not known, an in- 
direct relation connects the QCD picture with the hadron 
picture by QCD sum rules. In the hadron picture, the in- 
medium modifications arise from the coupling of mesons 
to resonance - hole states and the coupling of the modi- 
fied mesons to resonances. The best known example is the 
treatment of the A in the framework of the Z\-hole model 
(see e.g. Ref. [2J|3]). The hadron in-medium spectral func- 
tions for 7T-, 77-, and p-mesons and baryon resonances have 
been recently calculated by Post, Leupold, and Mosel [4] 
in a self-consistent coupled channel approach. 

The experimental investigation of hadron in-medium 
properties is complicated by initial and/or final state in- 
teractions. Since the present experiment uses photopro- 
duction of mesons, no initial but significant final state in- 
teraction effects must be considered. Here, the investiga- 
tion of these reactions also allows us to perform a detailed 
study of the meson - nucleus interactions which are re- 
sponsible for the final state interaction [SB?]. In case of 
the short-lived r\ meson the investigation of final state in- 
teraction effects is almost the only possibility to study the 
77-nucleon interaction. 

Experimentally, one of the clearest, although still not 
fully explained, in-medium effects has been observed in 
the excitation function of the total photoabsorption reac- 
tion [8]l9l[T0] . The bump in the elementary cross sections 
around 700 MeV incident photon energy, corresponding 
to the second resonance region, namely the excitation of 
the Pn(1440), Di 3 (1520), and Sn(1535) resonances, is not 
seen in the nuclear data. Many different effects have been 
discussed in the literature including trivial explanations 
like nuclear Fermi motion. Fermi motion certainly con- 
tributes to the broadening of the structure but cannot ex- 
plain its complete disappearance. Collisional broadening 
of the resonances due to additional decay channels like 
NN* — > NN has been studied in detail in the framework 
of transport models of the Boltzmann-Uehling-Uhlenbeck 
(BUU) type (see e.g. [11]) and can also not fully explain 
the data. The situation is complicated by the fact that al- 
ready on the free nucleon the second resonance bump con- 
sists of a superposition of reaction channels with different 
energy dependencies [6j . Inclusive reactions like total pho- 
toabsorption do not allow to study in-medium properties 
of individual nucleon resonances. A study of the partial re- 
action channels is desirable, but their experimental identi- 
fication is more involved, and final state interaction effects 
[7] as well as experimental bias due to the averaging over 
the nuclear density [12] must be accounted for (see Ref. 
|13] for a recent summary). Of special interest are meson 
production reactions which are dominated in the energy 
region of interest by one of the three resonances. Single 
and double pion production reactions have been employed 
for the study of the D13 resonance in the nuclear medium 
[51114). although up to now without conclusive results. 

Photoproduction of 77 mesons in the second resonance 
region is an excellent tool for the study of the Sn(1535) 
resonance, which completely dominates this reaction [15[ 
[TB] . Photoproduction of 77 mesons has been studied for the 
free proton in great detail over a wide range of incident 



photon energies and for different observables [Il) ] ll"7 ] [r8 . 19] 
[20 ] l2i p 22 p 23 j fe4]. The quasi-free reaction off the neutron 
bound in light nuclei has been investigated in detail for 
incident energies up to the peak position of the Su(1535) 
(« 800 MeV) [25,26,27], quasi- free neutron/proton cross 
section ratios for a few angular ranges up to photon en- 
ergies of 1 GeV have been reported in [28] and the co- 
herent photoproduction off light nuclei has been inves- 
tigated for the dcutcron and Helium isotopes 26,28,29, 
[30] . The combined result of these experiments (see [31] 
for a summary) was, that up to photon energies of «900 
MeV also on the neutron the reaction is completely dom- 
inated by the Sn(1535) with a constant cross section ra- 
tio cr n /<Jp « 2/3. Only very recently, results from the 
GRAAL, ELSA, and Tohoku experiments [32 ] [33 ]l34ll35] 
indicated a stronger contribution of a higher lying reso- 
nance to 771 — * 7771 than to 77J — > r/p for photon energies 
above 1 GeV. 

A first search for possible in-medium effects on the Sn 
spectral function was done with the TAPS experiment at 
MAMI |5j. However, the experiment covered only incident 
photon energies up to 800 MeV, i.e. approximately up to 
the peak position of the resonance. The experimental re- 
sults were in good agreement with BUU-model calcula- 
tions (see e.g. QJ]). Subsequently, measurements at KEK 
[36] and Tohoku [37] extended the energy range up to 
1.1 GeV. The KEK experiment reported some collisional 
broadening of the Sn resonance. The Tohoku experiment 
pointed to a significant contribution of a higher lying res- 
onance to the 771 — > 77,77 reaction. However, none of these 
experiments covered the full line shape of the Sn. 

Here, we report the measurement of 77 photoproduc- 
tion off carbon, calcium, niobium, and lead nuclei up to 
incident photon energies of 2.2 GeV, i.e. throughout and 
beyond the Sn resonance range. For comparison, the re- 
action has been studied for the same energy range off deu- 
terium which provides an estimate for the average nucleon 
cross section. 



2 Experimental setup 

The experiment was performed at the electron stretcher 
accelerator facility ELSA [38]l39] in Bonn, using a 2.8 GeV 
electron beam. Real photons were produced by Brems- 
strahlung off a copper foil of 0.3 % radiation length thick- 
ness. The photon energies were determined via the mo- 
mentum analysis of the scattered electrons by a mag- 
netic spectrometer. The tagging system, which is operated 
in coincidence with the production detector viewing the 
targets, is shown in Fig. [1] The direct electron beam is 
stopped in a beam dump while electrons having emitted 
Bremsstrahlung are deflected into the detection system of 
the tagging facility. The system has 14 overlapping scin- 
tillator bars with 4 cm thickness which cover the photon 
energy range between 22 % to 95 % of the incoming elec- 
tron beam energy E Q . Better energy resolution is provided 
by a scintillating fiber detector which covers 18 % to 80 % 
of E Q and a wire chamber (80 % to 92 %). In the present 
experiment only the scintillating fiber detector was used 
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Fig. 1. Setup of the tagging spectrometer. 




which provides an energy bin width of «10 MeV for the 
lowest incident photon energies around 650 MeV and 2 
MeV at the high energy end of 2.2 GeV. The total rate in 
the tagging system was 8-10 MHz for an incident electron 
beam intensity of ~ InA. 

Solid targets of 12 C (20 mm length), 40 Ca (10 mm), 
93 Nb (1 mm), and ™ Qt Pb (0.64 mm) were irradiated by the 
photon beam. The lengths of the carbon, calcium, and lead 
targets corresponded to 8 - 10 % of the respective radia- 
tion length Xq. The niobium target was somewhat thicker 
(« 17 % of Xq), all targets were 30 mm in diameter. The 
ry-mesons produced in the photonuclear reactions were de- 
tected via their rj — > 3tt° — > 67 decay (branching ratio 
32.5 %) with a two-component electromagnetic calorime- 
ter, covering 99 % of the full solid angle (see fig. [2]). The 
targets were mounted in the center of the Crystal Barrel 
detector [40] which covered the full azimuthal angle for 
polar angles between 30° and 168°. The Barrel consisted 
of 1290 CsI(Tl) crystals of 16 radiation lengths Xq. Inside 
it, around the target, a three-layer scintillating fiber detec- 
tor [41] (513 fibers of 2 mm diameter, three layers oriented 
with respect to the z-axis by angles of -24.5°, +25.7°, 
0°) was mounted for charged particle identification. Com- 
pared to the standard setup of the Crystal Barrel which 
was used for the measurement of ^-photoproduction off 
the proton [21] (see [42] for a detailed description of the 
setup), the 90 forward- most crystals have been removed. 
The forward angular range down to 4.5° was covered by 
the TAPS detector [431144] . This component consisted of 
528 BaF2 crystals of hexagonal shape with an inner diam- 
eter of 5.9 cm and a length of 25 cm corresponding to 12 
radiation lengths. They were arranged in a wall-like struc- 
ture as shown in the lower part of fig.[H A 5 mm thick plas- 
tic scintillator was mounted in front of each BaF2 crystal 
for the identification of charged particles. The front face 
of the BaF2 wall was located 1.18 m from the center of 
the target. Both calorimeters have a comparable energy 
resolution of [4"UIl4"4] 
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The impact points of photons are determined from the 
center of gravity of the electromagnetic showers, so that 



Fig. 2. Arrangement of the Crystal Barrel and TAPS detec- 
tors. Upper part: side view, lower part: front view of the TAPS 
wall: left hand side: logical segmentation for the LED-low trig- 
ger, right hand side: logical segmentation for the LED-high 
trigger (see text). 



the angular resolution is better than the granularity of 
the crystals. It is 1.5° (a) for the CB [40] for photons 
with energies above 50 MeV and 1.25° in TAPS. The fast 
BaF2 modules were read out by photomultipliers, the Csl 
crystals by photodiodes. Therefore, only information from 
the TAPS wall could be used for the first level trigger. For 
this purpose each module of the TAPS wall was equipped 
with two independent leading edge discriminators which 
were combined in two different ways into logical groups 
(see Fig. [2]). For the present experiment the thresholds of 
the first set of leading edge discriminators were set to 60 
MeV (LED-low) and the thresholds of the second set to 80 
MeV (LED-high). A valid first level trigger was accepted 
if either at least two logical groups of the low-threshold 
or at least one group of the high-threshold discriminators 
had fired. In the latter second level trigger from 

the FAst Cluster Encoder (FACE) of the Crystal Barrel, 
indicating at least two separated hits in the Crystal Barrel, 
was required in addition. Due to the trigger conditions 
only the decay channel into six photons could be used for 
the detection of the 77 mesons since the probability to find 
both photons from a two photon decay in TAPS is almost 
negligible. It should be noted that this restriction occurs 
only for measurements off nuclei where the recoil nucleon 
can be a neutron. In case of a proton target the recoil 
proton can provide the trigger. 
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3 Data analysis 

In the experiment, 77-mesons have been identified via the 
decay chain r\ — > 7r 7r°7r° — > 67. Events with six detected 
photons without a condition on further detected charged 
and/or neutral particles (recoil nucleons, pions) were se- 
lected. The photon reconstruction and identification in 
the Crystal Barrel is discussed in detail in Ref. [42]. It 
is based on a cluster search algorithm and uses the infor- 
mation from the three layer scintillating fiber detector for 
rejection of charged particles. The photon identification in 
TAPS is based on the information from the charged par- 
ticle veto detectors, on a time-of-flight analysis, and on a 
pulse-shape analysis of the BaF2 signals. Details of this 
analysis procedure can be found in Ref. 45j. 
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Fig. 3. Upper part: Invariant mass of photon pairs from events 
with six photons. All possible disjunct combinations of photons 
are included. Solid line: data, dashed line: Monte Carlo simu- 
lation of r\ — > 37T° — > 67. Bottom part: same data samples but 
only 'best' combination selected by x 2 test is plotted (see text). 
The vertical lines indicate the cut applied to select candidates 
for the 77 — > Sir decay 
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Fig. 4. Identification of ?y-mesons via invariant mass. Upper 
row: incident photon energies in the range 0.8 - 0.85 GeV, mid- 
dle: 1.4 - 1.6 GeV, bottom: 1.85 - 2.15 GeV. Left column: spec- 
tra with fitted background and rj-peak. Right column: back- 
ground subtracted peaks compared to Monte Carlo line-shapes. 
Vertical lines: nominal position of r] mass peak. All events have 
been corrected for detection efficiency (see text). 



events with six (or more) photons the invariant mass of 
all possible disjunct photon pairs was calculated. In the 
case of six photons 15 different combinations into pho- 
ton pairs are possible. The resulting two-photon invariant 
mass spectrum is shown in Fig. 2] (top part) together with 
a Monte Carlo simulation of the r\ — » 3tt° — > 67 decay. 
Among the 15 combinations the 'best' combination was 
chosen via a \ 2 analysis which minimizes the following 
expression 



Even in the presence of intense charged particle back- 
ground, the analysis leads to a very clean photon sample. 
The following invariant mass analysis is based on the ex- 
cellent reproduction of the shapes of the invariant mass 
peaks by Monte Carlo simulations. In the first step for 



E 

k=l 



(ro 77 fa, j k ) - m°) : 
(Am^(i k , j k )) 2 



(2) 



where ii, ...,13, ji, ...j'3 represents a permutation of 1,...,6, 
is the pion mass, and the the invariant masses 

of the photon pairs with their uncertainties Z\m 77 . The 
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resulting spectra for data and simulation are shown in 
Fig. 2] (bottom part). In case of the simulation, where 
the background in the upper part of the figure is only 
of combinatorial nature, a clean, background-free 7r° in- 
variant mass peak is recovered by this procedure. For the 
data, some background from other reactions remains. Due 
to the selection procedure of the 'best' combination this 
background is also concentrated around the peak region. 
In the next step, events were selected where all three two- 
photon invariant masses of the best combination are lying 
between 110 MeV and 160 MeV (indicated by the ver- 
tical lines in fig. [4|. This cut is motivated by the shape 
of the simulated invariant mass cut and removes only a 
small fraction of 'true' events, which is determined from 
the simulation and taken into account for the extraction 
of the cross section. 

The nominal mass of the pion was then used as a con- 
straint to improve the experimental resolution. Since the 
angular resolution of the detector for photons is much bet- 
ter than the energy resolution, it was not necessary to use 
a kinematic fitting procedure. Instead, only the photon 
energies were re-calculated from 



E'i.2 — El, 



(3) 



where E\_2 ar e the measured photon energies, E[ 2 the 
re-calculated, m T o is the nominal ir° mass, and m 77 the 
measured invariant mass. This procedure was applied to 
all three photon pairs combining to pions. In the last step, 
the invariant mass spectrum of the six photons is build 
from their four-vectors, using the measured angles and the 
re-calculated energies. Typical results for the rj invariant 
mass peak for different incident photon energies are show 
in Fig. 21 The spectra can be fitted by a simple polyno- 
mial background (polynomial of second degree) and the 
line shape of the invariant mass peaks generated from a 
Monte Carlo simulation with the GEANT3 package [46] . 
For the fit only the three parameters of the background 
polynomial and the amplitude of the simulated response 
function were varied. Background subtracted spectra com- 
pared to the simulated line-shape are also shown in Fig. 
[4j The shape of the invariant mass peaks is in excellent 
agreement with the results of the Monte Carlo simulation 
(see Fig. where the same analysis procedure was ap- 
plied. The position of the invariant mass peak as function 
of incident photon energy is rather stable, which is partly 
due to the fact that the positions of the ir° and rj invari- 
ant mass peaks have been used in an iterative procedure 
for the energy calibration of the detector, and partly due 
to the re-calculation of the photon energies from the ir° 
invariant masses. For differential cross sections, this pro- 
cedure must be applied to each bin. The detection effi- 
ciency, as discussed in the next section, was corrected on 
an event-by-event basis. Therefore, the fitting was not ac- 
tually done on the raw invariant mass distributions but on 
the corresponding spectra with efficiency corrected events 
(shown in fig. |4|). 

With the analysis discussed so far, the fully inclusive 
reaction jA — > rjX is identified, where X may also in- 




200 400 

Am [MeV] 

Fig. 5. Missing mass spectra for lead calculated under the 
assumption of quasi-free single r\ production. Left hand side: 
incident photon energies in the range 0.65 - 0.8 GeV, right hand 
side: 1.3 - 1.5 GeV. Histograms: simulated detector response 
for quasi-free single rj production. Arrow indicates cut for this 
reaction. 



elude any (kinematically possible) number of pions. Se- 
lecting exclusive, single reproduction without any further 
mesons in the final state is not possible by simply vetoing 
events with additional clusters in the detector. Additional 
charged mesons may go undetected (e.g. along the beam 
line), making the suppression incomplete. In case of the 
Barrel, charged pions cannot be distinguished from pro- 
tons and the condition would falsely suppress events with 
detected recoil protons. Therefore, exclusive //-production 
can only be identified via the reaction kinematics. Here, it 
is assumed that the reaction occurs quasi-free off a bound 
nucleon. The initial momentum of the nucleon is neglected 
and the missing mass Am of the reaction is calculated 
from the nucleon mass tun, the energy of the incident 
photon Eb, and the energies and momenta E li and P~ ii of 
the six decay photons: 



Am = 



\ 



(E b + m N 



6 



m N . 



i=l 



(4) 

The resulting distributions are broadened by Fermi mo- 
tion, so that a perfect separation of the different reac- 
tion channels is not possible. Examples of missing mass 
spectra are shown in fig. [5] The structures around zero 
missing mass are related to single ^-production, the con- 
tribution at large positive missing mass, which is only vis- 
ible at higher incident photon energies, originates from 
rjir final states and secondary production processes like 
"/N — » 7riV, nN — > rjN. The experimental results are com- 
pared to a simulation of the detector response based on the 
shape of the missing mass distributions for single quasi- 
free ^-production predicted by a BUU model (see sec.[H]). 
The indicated cut was used for the analysis of single re- 
production. This cut does not allow a perfect separation of 
single ^-production from the other contributions since the 
tails of the different distributions are overlapping. How- 
ever, since the same cut was used for the model results 
(see sec. H]), the comparison between data and model is 
not affected. 
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4 Determination of cross sections 

The absolute normalization of the measured yields was ob- 
tained from the target densities, the incident photon flux, 
the 77 — * 3ir° — ► 67 decay branching ratio (b^e^— 31.35 
%), and the detection efficiency of the calorimeter. The 
measurement of the photon flux is based on the count- 
ing of the deflected electrons in the focal plane detectors 
of the tagging spectrometer (see. fig. []}. The fraction of 
correlated photons, passing the collimator and impinging 
on the target, was determined roughly once per day in a 
mode where the trigger is derived from the tagger and the 
photon flux (at reduced beam intensity) is measured by 
a photon counter placed downstream of the calorimeter. 
The detector dead-time of approximately 60 % was deter- 
mined with scalers gated by life-time of the experiment 
and spill-time of the accelerator, respectively. 




does not vanish for any kinematically possible combina- 
tion of 77 polar angle and kinetic energy. This is demon- 
strated in fig. [6] where the correlation between angles and 
energies of the measured 77 mesons is compared to the 
simulated detection efficiency. The absolute values of the 
detection efficiency are not large since the first level trig- 
ger was only sensitive to photons in the TAPS forward 
wall (see sec. [5]). Even for the six-photon decay of the 77 
the efficiency of this trigger condition is small, in partic- 
ular for 77 mesons emitted at large polar angles. However, 
it is obvious from the figure, that the entire phase space 
of kinematically possible combinations is covered by non- 
vanishing detection efficiency, so that for the determina- 
tion of cross sections no extrapolations had to be done. 



5 Systematic uncertainties and comparison to 
previous results 

The main sources for systematic uncertainties are related 
to the background level in the ^-invariant mass spectra, 
the simulation of the detection efficiency, and the deter- 
mination of the incident photon flux. Other uncertainties 
like e.g. the surface thickness of the solid targets (better 
than 1 %) are comparably negligible. 



Fig. 6. Left hand side: laboratory angle and kinetic energy 
distribution of measured 77-mesons. Right hand side: detection 
efficiency as function of the same parameters (regions clearly 
outside kinematically possible combinations have not been sim- 
ulated) . 



The detection efficiency of the Crystal Barrel/TAPS 
calorimeter has been determined by Monte Carlo simu- 
lations using the GEANT3 package [46]. The simulation 
includes all relevant properties of the experimental setup 
like geometrical acceptance, trigger efficiency, detection 
efficiency, and analysis cuts. The branching ratio for the 
r/ — > 67 decay is not included in the efficiency. The 77 
mesons are produced in many different final states and 
reaction types: single 77, rjir, re-scattered 77 mesons, and 77 
mesons from secondary reactions in the nucleus. All pro- 
cesses are additionally complicated by the momentum dis- 
tribution of the bound nucleons. Consequently, the cor- 
relation between momentum and emission angle of the 
mesons is a priori not known. Therefore, a reliable model 
is not available for an event generator for the Monte Carlo 
simulations. Instead, the detection efficiency has been de- 
termined from the simulations as a function of the lab- 
oratory polar angle and laboratory kinetic energy of the 
77-mesons, which are measured quantities. It was then ap- 
plied on an event-by-event basis to the data. The method 
is described in more detail in [6] for inclusive 7r° produc- 
tion off nuclei. In this way, a model independent detection 
efficiency correction is achieved as long as the efficiency 
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Fig. 7. Typical systematic effects (carbon data). Main plot: 
total cross section from detection efficiency with 10° binning 
in 9 V (filled (red) circles) compared to 2° binning (open (black) 
circles). Bottom insert: ratio of cross sections. Top insert: typ- 
ical systematic uncertainty related to background shape in 77- 
invariant mass spectra (see text). 

As discussed in sec. |3l the background level beneath 
the 77 signal has been determind by fitting the amplitude 
of the simulated shape of the invariant mass peak and 
the three parameters of a background polynomial. The 
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Fig. 8. Comparison with earlier results. Upper part: 12 C from 
Mainz [5 , Tohoku [37] and this work. Bottom part: 40 Ca from 
Mainz [45], 63 Cu from Tohoku [37], and 40 Ca and 93 Nb this 
work. The error bars of the present data include only the statis- 
tical errors, the shaded bands at the bottom indicate the total 
systematic uncertainties for 12 C (top) and 40 Ca (bottom). 



resulting background subtracted signal, shown in fig [4] 
agrees well with the simulated line shape. The systematic 
uncertainty of this procedure was studied by a variation 
of the fitted range in the spectra, giving rise to systematic 
deviations in the background shape. Typical systematic 
variations of the fitted peak amplitude (see fig. [7] upper 
insert) range from 2% at low incident photon energies to 
8% at the highest incident photon energies. 

Since the detection efficiency could be simulated with- 
out any assumptions about the angular and energy dis- 
tributions of the mesons, the corresponding uncertainties 
are small, estimated at the 5 % level. They are mainly due 



to the exact representation of thresholds, shower develop- 
ment, absorbing inactive material in the detector, and the 
exact target and beam positions in the GEANT simula- 
tions. The stability of the detection efficiency correction 
has been checked by a variation of the bin size of the 
simulated efficiency by a factor of five (Fig. [6] shows the 
coarsest binning in angle) . Figure [7] shows a comparison 
of the total cross section for carbon constructed with bin 
sizes of 10° and 2°. Most data points from the two anal- 
yses agree within ±2.5% and the fluctuations are mostly 
within statistical uncertainties (see fig. [7] bottom insert). 

The photon flux involves the determination of the rate 
of scattered electrons and the measurement of the tagging 
efficiency. The uncertainty could be estimated from the 
measurement with the deuteron target. In this case, dif- 
ferent incident photon spectra have been used. About half 
of the data were measured with an electron beam energy of 
2.6 GeV and a linearly polarized photon beam (coherent 
Bremsstrahlung from a diamond lattice, see |24j for de- 
tails) with the polarization maximum around 1 GeV. The 
other data were measured with an electron beam energy 
of 3.2 GeV and unpolarized photon beam. Consequently, 
the shapes of the incident photon spectra were different, 
due to the coherent peak. Due to the different beam ener- 
gies, scattered electrons corresponding to the same photon 
energy have been registered in different sections of the tag- 
ging spectrometer. After the flux correction, the cross sec- 
tions for r\ production extracted from the two runs agreed 
within ±10 %, taken as the typical systematic uncertainty. 

The flux uncertainty is identical for all nuclei (same 
settings of incident photon energy, beam intensity, and 
parameters of the tagging system for all nuclei) and thus 
does not influence A-dependent properties. It must only be 
accounted for in the comparison to model results or data 
from other experiments. For this case, all three systematic 
effects have been added in quadrature. The systematic un- 
certainties of the detection efficiency and the background 
subtraction are also not completely independent for the 
four nuclei. However, for an estimate of systematic uncer- 
tainties of scaling properties we have made the most pes- 
simistic assumption that this effects vary independently. 

The total cross section data for inclusive rj production 
are compared in fig. [8] to previous results below 1.2 GeV. 
Carbon and Calcium data are available from Mainz below 
0.8 GeV 0[45], Carbon data from KEK [36] and Carbon 
and Copper data from Tohoku [37] below 1.2 GeV. In 
the Carbon case, the KEK and Tohoku data are system- 
atically higher by roughly 10 %. For the heavier nuclei, 
a direct comparison between the Tohoku results and the 
present data is not possible, since different target nuclei 
have been investigated. At incident photon energies below 
800 MeV, all data scale like A 2 / 3 . At higher incident pho- 
ton energies, the A 2 / 3 scaling does not hold anymore, the 
present Calcium data are clearly lower than the present 
Niobium data (see also fig. [5]). From this behavior, one 
expects that the scaled cross section for Copper lies in 
between the Calcium and Niobium results. However, the 
Tohoku Copper results fall on top of our Niobium data, 
indicating that also for the heavier nuclei the Tohoku re- 
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suits are systematically higher by roughly 10 %. However, 
the energy dependence of the excitation functions for the 
present data and the KEK and Tohoku results is very 
similar. Rather good agreement is found when all data 
are re-scaled by 10 %, which is within their systematic 
uncertainty. 



6 The BUU-model calculations 

The detailed interpretation of the experimental results 
is only possible via a comparison to model calculations 
which incorporate effects like nuclear Fermi motion, Pauli- 
blocking of final states, and in particular the propagation 
and absorption of mesons and nucleon resonances in nu- 
clear matter. Results obtained in the framework of the 
BUU transport model for photon induced reactions as dis- 
cussed in detail in [4TPMTT] have been used. The model 
is based on the BUU equation: 

^ + V P H ■ V r - V r H ■ V p ^j Fi{r,p, t) = 

= Icou [Fn, F^, F v , Fn* , Fa, ...] (5) 

which describes the space-time evolution (r: space coor- 
dinate, p: momentum) of the spectral phases space density 
Fi of an ensemble of interacting particles i=N, N*, A, 7T, ?y, . 
in nuclear matter from the moment of their creation to 
their absorption or escape through the outer boundaries 
of the nucleus. The left hand side of the equation - the 
Vlasov term - describes the propagation of the particles 
under the influence of a Hamilton function H . It contains 
information about energy, mass, self-energy (mean field) 
of the particle and a term that drives back an off-shell 
particle to its mass shell. Explicitly, it can be written as 

H = v/Ou + S) 2 + p 2 (6) 

with the particle mass \i and a scalar potential S for 
baryons [IT] . The right hand side of the BUU equation 
- the collision integral - describes particle production and 
absorption. It consists of a gain and a loss term for the 
phase space density Fi, accounting for interactions be- 
tween the particles beyond the mean-field potential. 

The constituents of the nucleus are defined as 'test 
nucleons' and follow a Woods-Saxon density distribution 

where the nuclear radius is related to the nucleus mass 
via R = 1.124A 1 / 3 fm and a = (0.0244 A 1 / 3 + 0.2864)/m. 
The momentum distribution is described within the Fermi 
gas approach. 

/o 2 \V3 

Mr) = [— P(r)) ■ ( 8 ) 

The elementary r\ cross sections off protons and neu- 
trons are included in this model. The produced resonances 
and mesons propagate in the nucleus and can be scattered, 



absorbed or decay. The different reaction probabilities are 
either fitted to experimental data or calculated. They are 
incorporated into the model by the collision term and may 
interact according to the geometrical condition that the 
distance between the two particles is smaller than the im- 
pact parameter b c = y 'a /tt where a is the reaction cross 
section. 



7 Results 

The total inclusive r\ production cross sections (i.a. for 
the reaction — > r]X, without any condition for X) are 
summarized in fig. [9] At incident photon energies below 
«800 MeV the cross sections scale with A 2 / 3, (A atomic 
mass number) for the heavier nuclei and agree with the av- 
erage nucleon cross section (deuteron cross section scaled 
by a factor of 2). In the following A e ft means A=2 for 
the deuteron and A 2 ! 3 for all other nuclei. This behavior, 
which indicates strong absorption of the mesons, was al- 
ready found in [6] for rj, 7r°, and double tt photoproduction 
in the same energy range. However, at higher energies, the 
cross sections behave completely differently and scale al- 
most with the mass number. 

This is shown more quantitatively in fig. [9] (bottom) 
where the scaling coefficient a obtained from a fit of 

a(A) oc A a (9) 

is plotted versus the incident photon energy. It would be 
tempting to argue that the rise of the scaling coefficient 
simply reflects a decrease of the absorption cross section 
with increasing kinetic 77-energy, since the most efficient 
absorption process is s-wave excitation of the Sn(1535) 
resonance. However, the situation is not that simple. For 
the further discussion, we must keep in mind, that the 
scaling is not only influenced by the absorption cross sec- 
tion of the ^-mesons, but may also reflect A-dependent 
differences of their production. In the most simple case of 
quasi-free single-meson production, the production before 
FSI will scale with A and then a scaling of the observed 
meson rates with a=2/3 indicates strong FSI, while a scal- 
ing with A indicates transparent nuclear matter. In this 
case, the scaling coefficient will only depend on the ki- 
netic energy T of the mesons (energy dependent absorp- 
tion cross section). 

However, the data show a different behavior. As shown 
in Fig. [TU] for fixed values of T the scaling is dependent 
on the incident photon energy E 1 . Furthermore for fixed 
incident photon energy the coefficients drop from values 
close to unity for small T to roughly 2/3 for the largest 
T possible at that incident photon energy. This is exactly 
the opposite of what one would expect when the behav- 
ior of the scaling coefficients would be dominated by the 
s-wave absorption into the Sn. However, the observed be- 
havior can arise when the production rates of the mesons 
before absorption do not scale with A. Problematic in this 
respect are side-feeding contributions from secondary pro- 
duction processes like 7A — > ttN, ttN — > rjN, and FSI 
processes that modify the observed energy distribution of 
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Fig. 9. Upper part: Total fully inclusive rj cross section nor- 
malized by A eff = A 2/3 for A = 12,40,93,208 and A eff = 2 
for the deuteron. Error bars are statistical, the shaded bands 
at the bottom show the systematic uncertainties for carbon 
(lower double shaded band) and lead (upper band) excluding 
the 10% flux normalization uncertainty. Bottom part: Scaling 
coefficient a (see text) as function of incident photon energy. 
Black circles: results if only statistical uncertainties are con- 
sidered. Slightly displaced open symbols: including systematic 
uncertainties except 10% overall normalization. Curve: result 
from BUU model. The two inserts show the scaling with mass 
number for two typical ranges of incident photon energy. 



the rj mesons. They may completely obscure the effects 
related to ^-absorption, since they may strongly increase 
with mass number. The BUU-modcl simulations predict, 
that both contributions are substantial (see sec. [5]). 

However, due to energy and momentum conservation, 
secondary production processes (as well as rjir final states) 
cannot contribute at the kinematical limit (i.e. at maxi- 
mum T for given Ey), but will produce rj mesons with 
smaller kinetic energies (some energy is carried away by 
the additional nucleon(s) involved in the secondary reac- 



tions). Therefore, it is possible to extract the ?7./V-absorption 
cross section from the scaling behavior in this regime. For 
this purpose, the scaling factors a have been fitted for 
the high energy end of the T distributions for different 
incident photon energies using the condition 



T > (E 7 - m„)/2 



(10) 



where E 7 is the incident photon energy and the mass 
of the rj mesons. The result as function of rj kinetic energy 
is compared in fig.[Tl]to the previous results for low-energy 
•q mesons [5] and for 7r° mesons [6] . The somewhat surpris- 
ing result is that for rj mesons the scaling coefficient a is 
almost constant at 2/3 over a large range of rj kinetic en- 
ergy, indicating strong absorption independent of kinetic 
energy. In the case of pions, the absorption is expectedly 
very weak for kinetic energies, which are too low to ex- 
cite the A resonance. The pions escape from the nucleus 
and a (pions) is one. The scaling factor reaches 2/3 in 
the A regime and then seems to slowly increase again. 
The large absorption cross section for rj mesons at small 
kinetic energies was expected since this is the excitation 
region of the Sn resonance with a strong coupling to the 
Nrj channel. Unexpectedly, a decrease of the absorption 
cross section is not observed, even at kinetic energies far 
above this range. The corresponding rjN absorption cross 
section can be deduced from the results of a Glauber 
model calculation discussed in [5]. The model is based on 
the assumption, that secondary production processes play 
no role, which has been assured, as discussed above, by the 
choice of the kinematical conditions. The dependence of 
the scaling coefficient a on is shown in fig.H21 It yields 




500 750 

T[MeV] 

Fig. 10. Evolution of the scaling factor a with the kinetic 
energy for different ranges of incident photon energies. The 
open, slightly displaced symbols show for the highest incident 
photon energy the result when systematic uncertainties except 
the 10% flux normalization are included. 
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Fig. 12. Dependence of the scaling coefficient a on the rjN 
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a nN ~ 3® mD for the average value of a « 0.66. In a sim- 
ilar analysis, recently Kotulla et al. [49] have investigated 
the scaling behavior of the photoproduction of w-mesons 
off nuclei. They found typical absorption cross sections in 
the range of 50 mb, corresponding to inelastic in-medium 
widths of the mesons around 150 - 200 MeV. Here, we do 
not folluw-up this analysis quantitatively, however, it is 
evident that also the extracted 77 absorption cross section 
must correspond to inelastic widths at least in the few ten 
MeV range. 

Due to the contribution from tjtt final states and sec- 
ondary production processes, the inclusive reaction can- 
not be used to extract an in-medium line shape of the 
Sn(1535) resonance. A separation of quasi-free single 77 
production can only be achieved by cuts on the reac- 
tion kinematics. A cut on the missing mass (see fig. [5]) 
at Am < 140 MeV is motivated by a comparison of the 
missing mass spectra to the simulated line shape of the re- 
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Fig. 13. Total r\ cross section with missing mass cut at 140 
MeV. Errors are statistical, the shaded band shows the sys- 
tematic uncertainty (excluding the 10% flux normalization) for 
Calcium. 



sponse for quasi-free single r\ production. The total cross 
sections after the cut are summarized in fig. 1131 

The shape (position and width) of the Sn-resonance 
structure is very similar for carbon, calcium, niobium, and 
lead. A clear systematic evolution with mass number is not 
observed. The shape is different for the deuteron data but 
this effect can mostly be explained by the different mo- 
mentum distributions of nucleons bound in the deuteron 
or a heavy nucleus. 

With the exception of the deuteron target, the separa- 
tion of single quasi-free r\ production from other processes 
with the missing mass cut is only an approximation, due 
to the overlapping tails of the distributions from different 
processes. However, it gives already an indication, that 
no strong in-medium effects on the shape of the Sn res- 
onance occur. A more detailed discussion is possible by a 
comparison to the results of the BUU-model subjected to 
the same kinematical cuts. 



8 Comparison to BUU results 

The distributions of kinetic energy and cm-angle (cm sys- 
tem of the incident photon and a nucleon at rest) for the 
inclusive data are compared in Fig. [T3] to the results of 
the BUU model. The overall agreement between data and 
model is quite good. The most significant disagreement is 
observed for the differential cross sections at small kinetic 
energies of the 77-mesons. In this regime, the model cal- 
culations are closer to the simple A 2 / 3 scaling and under- 
estimate the observed cross sections for the heavy nuclei. 
As discussed below, these discrepancy can be traced back 
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Fig. 14. Upper part: energy distributions for different inci- 
dent photon beam energy ranges. Bottom part: angular dis- 
tributions. Curves are BUU model calculations. Uncertainties 
include systematic effects except the 10% flux normalization. 



to the contribution from rjir final states and/or secondary 
production processes of 77-mesons. At low incident photon 
energies the angular distributions are similar to the mea- 
sured deuteron distributions (i.e. to the average nucleon 
coss section) and agree quite well with the model results. 
At the highest incident photon energies, the angular dis- 
tribution for the deuteron peaks at forward angles, since 
for the free nucleon t-channel processes become impor- 
tant. This effect is not seen for the heaviest nuclei, where 
the distributions peak at backward angles. The model re- 
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Fig. 15. Upper part: Comparison of the total exclusive sin- 
gle 77 production cross (missing mass cut at 140 MeV) to 
BUU results. Shaded bands: total systematic uncertainty. Bot- 
tom part: filled symbols: inclusive cross sections (solid curves: 
BUU results). Open symbols: difference of inclusive and exclu- 
sive single 77 production (dashed curves: BUU-results). Shaded 
bands: total systematic uncertainty of the open symbols. 



suits show the same tendency and in the model the back- 
ward peaking contribution arises mainly from secondary 
production processes. This behavior is easily understood 
since 77-mesons from secondary production processes on 
average have small kinetic energies and therefore appear 
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at backward angles in the fast forward moving photon - 
nucleon cm-system. 

The total cross sections for the inclusive reaction, for 
quasi-free single 77 production and for the contribution 
from t/7t final states and secondary production processes 
are summarized and compared to the model results in Fig. 
[131 The shape of the total inclusive cross section is reason- 
ably well reproduced for the lighter nuclei, but disagrees 
significantly for lead, where the model shows still the peak 
of the Sn resonance, which is absent in the data. How- 
ever, this systematic shape change from light to heavy 
nuclei is not related to an in-medium modification of the 
Sn resonance. This is clearly demonstrated by the sepa- 
ration of the inclusive cross section into quasi-free single 
^-production and the other components. The separation 
has been done by the missing mass cut at Am < 140 MeV, 
which has been applied to data and model calculations (to 
the latter after folding them with the experimental resolu- 
tion). Although, as discussed above, this separation is not 
perfect, the result for the line-shape of the Sn dominat- 
ing quasi-free single 77 production is clear. Position, width, 
and peak cross section of the S11 agree for all nuclei quite 
well with the model results. Only the peak cross section 
shows a little systematic evolution from carbon (slightly 
overestimated) to lead (slightly underestimated), which is 
however within the systematic uncertainty introduced by 
the missing mass cut. The shape is fully explained by the 
'trivial' in-medium effects included in the BUU model, 
in particular the momentum distribution of the nucleons 
bound in heavy nuclei. 

The agreement is less good in the energy range above 
the Sn resonance, where the model overestimates the mea- 
sured cross sections. However, one must keep in mind, 
that here the separation by the missing mass cut becomes 
strongly dependent on the exact shape of the missing 
mass distributions for the different components. Indeed 
the main effect at higher incident photon energies seems 
to be a strong underestimation of the contribution from rjir 
final states and/or secondary processes (see Fig. bot- 
tom part). In particular for lead around 1 GeV incident 
photon energy, this component rises much more rapidly in 
the data than in the model. 

This mismatch is most clearly seen in the missing mass 
spectra in this energy range which are compared in fig. 
[TBI to the BUU calculation. The contribution from single, 
quasi-free 77 production is overestimated while the contri- 
bution from ?77r final states seems to be underestimated in 
the model. Part of this discrepancy is probably due to the 
uncertainty in the elementary cross sections for rjir pro- 
duction reactions. There are recent precise data for the 
773 — > pr]TT reaction [50j . however, much less is known for 
the channels with charged pions in the final state or a 
neutron in the initial state. In this range of incident pho- 
ton energy, the modeled missing mass plots seem to indi- 
cate that the quasi-free single 77-peak has already signifi- 
cant contamination from the tails of the 'background' pro- 
cesses, in particular, from secondary reproduction. How- 
ever, it is not possible to obtain a reasonable fit of the mea- 
sured missing mass distributions by a variation of the ar- 
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Fig. 16. Missing mass spectra for carbon and lead for two 
ranges of incident photon energies. Curves: BUU-results for 
full model (full curves), single, quasi- free ^-production (dot- 
ted), 7/7r final states (dash-dotted) and secondary 77- production 
(dashed) . 



eas of the three model contributions, keeping their shape. 
Fitting the part of large missing mass with the r/ir and 
secondary 77 production contributions leads to unreason- 
able contributions of their tails in the quasi-free region 
around zero missing mass. Therefore not only the magni- 
tude but also the shape of this contributions seems to be 
partly in conflict with the data. On the other hand, the 
missing mass shape of the quasi-free single 77 production 
seems to be in better agreement with the data (see Fig. 
[SJ, it certainly agrees with it below the rjir threshold. 



9 Conclusions 

The investigation of inclusive and exclusive 77 production 
cross sections for heavy nuclei from threshold to 2 GeV 
can be summarized as follows. In the excitation region of 
the Sn (1535) resonance, contributions from ?77r final states 
and secondary production processes to inclusive 77 pro- 
duction are already significant. At higher energies these 
contributions even become dominant. A discussion of in- 
medium properties of the Sn resonance or of absorption 
properties of 77 mesons in nuclear matter requires a careful 
treatment of these effects. 

An analysis of the scaling of the cross sections with 
atomic mass number has been performed for rj mesons pro- 
duced closely to the kinematical limit where only quasi- 
free single 77 production can contribute. Combined with 
previous low energy results [5] , it is found that the scaling 
coefficient a is almost constant at a value of 2/3 for 77 ki- 
netic energies from 20 MeV up to 1 GeV. Using a simple 
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Glauber model approximation, this corresponds to a con- 
stant ijN absorption cross section of ~30 mb. A decrease 
of the absorption probability for ry-mesons with kinetic 
energies much above the Sn range is not observed. 

An analysis of the line shape of the Sn resonance can 
be achieved with the results for single, quasi-free r\ pro- 
duction after cuts on the reaction kinematics. The ob- 
served excitation functions for heavy nuclei have almost 
identical shape from carbon to lead. The results in the 
Sn range are in good agreement with BUU model calcu- 
lations which include the 'trivial' in-medium effects like 
Fermi smearing, Pauli blocking of final states, and contri- 
butions from secondary processes. Thus, an indication of 
a shift or broadening of the resonance has not been found. 

At higher incident photon energies, the agreement be- 
tween BUU calculations and experiment is less good. The 
relative contribution of single, undisturbed 77 photopro- 
duction is overestimated in the model and the contribu- 
tion of secondary processes and/or t/tt final states is sig- 
nificantly underestimated. This indicates a need for better 
input for the semi-inclusive rjX channels in the BUU cal- 
culations. 
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